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The  effects  of  solar  systems  operation  can  be  compared  only  under  very  similar  weather  conditions. 
Diagnostics  of  the  solar  systems  requires  unequivocal  determination  of  solar  irradiation.  Development 
of  a  method  for  precise  identification  of  solar  radiation  day  time  profile  is  needed,  as  the  methods  used 
so  far  in  the  cloud  cover  determination  are  not  satisfactory.  The  paper  presents  two  optional  methods, 
developed  by  the  authors,  for  identification  of  the  solar  radiation  profile.  Advantages  and  disadvantages 
of  the  methods  are  also  specified. 
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1.  Introduction 

Accurate  assessment  of  daily  operational  efficiency  of  the  solar 
thermal  heating  and  photovoltaic  (PV)  systems  calls  for  unequiv¬ 
ocal  determination  of  a  curve  that  characterize  distribution  of  the 
daily  sun  radiation.  This  curve  is  also  needed  for  comparative  anal¬ 
ysis  of  many  other  thermal  systems. 

Previous  studies  on  similar  problems  have  been  concentrated 
on  the  summary  daily  sun  radiation  and  sometimes  on  cloudiness 
degree  as  the  key  factors  to  determine  operational  efficiency  of  the 
examined  solar  systems.  Thermal  efficiency  is  assumed  as  directly 
proportional  to  daily  solar  irradiation,  on  condition  that  it  is  greater 
than  its  threshold  value  (e.g.  ISO  9459  [1]).  To  some  degree  it  is  a 
simplified  assumption.  Sun  radiation  intensity  forces  modification 
of  its  distribution  curve.  Irradiance  variability  in  time  determines 
the  heat  transfer  dynamics.  Every  change  in  distribution  curve  is 
critically  important,  as  it  greatly  affects  the  assessment  results  of 
solar  thermal  systems’  operational  efficiency. 

Radiation  variability  at  a  given  geographical  location  results 
from  two  overlapping  phenomena:  (1)  parabola  of  sun  movement 
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on  the  sky  is  determined  by  a  subsequent  day,  and  (2)  a  stochastic 
course  of  cloudiness.  If  sky  is  clear,  it  is  relatively  easy  to  determine 
variability  of  total  radiation  intensity  for  the  determined  atmo¬ 
sphere  turbidity  [2]. 

Many  empirical  studies  indicate  the  need  for  more  precise  iden¬ 
tification  of  the  radiation  course,  not  only  on  the  basis  of  average 
hourly  values  [3,4].  Assessment  of  daily  operational  efficiency  of 
solar  heating  and  photovoltaic  systems  should  be  based  on  data 
that  characterize  irradiance  changeability  in  the  short  intervals  (in 
hours  and  even  shorter  periods). 

The  authors  of  this  paper  propose  two  empirical  methods 
for  parameterization  of  the  sun  radiation  variability,  including 
description  of  the  methods  and  the  results  of  their  practical 
application. 

2.  Classification  of  changes  in  solar  irradiance  conditions 

Cloudiness  of  stochastic  character  determines  the  instantaneous 
irradiance  value.  In  meteorology,  the  cloudiness  is  determined  by 
firmament  coverage  with  clouds;  if  clouds  cover  the  sky  totally  the 
cloudiness  is  designated  with  8,  if  the  sky  is  totally  clear  then  0  is 
assigned.  Number  9  is  additionally  reserved  for  the  case  when  the 
cloud  cover  cannot  be  determined  because  of  smog,  fog,  etc. 
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The  formulation  of  solar  irradiance  models  is  a  complex  and 
difficult  task,  particularly  when  obtaining  of  needed  data  of  good 
quality  is  impossible.  In  such  cases,  the  cloud  coverage  observations 
are  used  as  basis  information  [5]. 

Belcher  and  DeGaetano  [6]  provided  information  on  typical 
hourly  cloud  cover  classes  for  the  United  States  conditions  and  the 
use  of  Automated  Surface  Observation  System  (ASOS).  They  for¬ 
mulated  the  solar  irradiance  model  and  determined  seven  types  of 
hourly  cloudiness: 


station  data  in  Catalonia  (Spain)  divided  every  hour  into  twelve 
5-min  intervals  N=  12,  where  average  solar  irradiation  value  was 
recorded  [15].  They  proposed  three  indexes  of  hourly  variability  for 
the  standardized  clearness  index  K'T  (standardization  is  connected 
to  current  value  of  optical  mass).  The  first  one  is  based  on  standard 
deviation  of  mean  value  KrT  of  the  whole  hour: 


Kf 


•  clr  -  clear  sky; 

•  few  -  few  clouds  on  the  sky; 

•  set  -  scattered  clouds; 

•  bkn  -  broken  clouds  (there  are  some  clear  sky  breaks  and  big 
cloud  groups); 

•  ovc  -  overcast; 

•  obs  -  total  cloud  cover; 

•  multi  -  hours  of  observations  when  at  least  one  cloud  has 
appeared. 


The  visual  evaluation  of  cloud  cover  made  by  individual 
observers  differed  distinctly  from  ASOS  results,  particularly  in  clr 
category  (almost  twice  more  frequent  hours  were  indicated  by 
ASOS)  and  in  the  set  and  bkn  categories  (over  twice  bigger  fre¬ 
quency  of  such  hours  were  indicated  by  visual  observers). 

The  cloud  cover  degree  and  complementary  clearness  index  I<T 
(the  ratio  of  global  horizontal  radiation  G  to  horizontal  extraterres¬ 
trial  radiation  G0)  and  sunshine  hours  are  basic  parameters  used  to 
formulate  solar  radiation  models  for  different  regions  of  the  world. 
Tens  of  correlation  models  have  been  formulated  all  around  the 
world,  to  examine  the  relationship  between  the  ratio  of  mean  daily 
global  irradiation  in  a  given  month  ( Gday )  and  monthly  average 
extraterrestrial  radiation  values  ( Go  day )  and  the  ratio  of  average 
daily  hours  of  bright  sunshine  ( Sday )  and  the  average  solar  day 
length  (Saday)  [7]: 


a  + 


^  x  $day 
$o,day 


The  numerical  values  of  the  regression  coefficients  a  and  b  differ 
when  other  independent  variables  including  location  latitude,  sun 
declination,  atmosphere  temperature,  etc.,  are  used. 

Development  of  measurement  technique  enabled  to  determine 
(characteristic  for  a  given  region)  the  classical  correlation  between 
the  share  of  diffuse  radiation  Dhor  in  total  radiation  sum  Ghor  (diffuse 
ratio)  on  a  horizontal  surface  and  the  clearness  index  Kj  [8]: 


Dhor 

Ghor 


( 


where  standard  deviation  has  the  form: 

a  = - fTi - 

The  second  index  contains  the  notion  of  absolute  deviation  sum 
between  subsequent  5-min  values: 

N-l 
)KT  l=1 

The  third  index  takes  into  account  the  difference  between  extreme 
values  in  a  given  hour: 

A3  =  ln  ((N-l)K^m“  '  ^min') 

Indexes  of  hourly  variability  of  clearness  index  and  its  mean  value 
were  then  used  (in  various  combinations)  in  multi-parameter  mod¬ 
els  to  describe  the  diffuse  irradiation  fraction.  The  results  in  the 
form  of  linear  least-squares  regressions  were  compared  to  analog¬ 
ical  dependences  obtained  according  to  the  classical  procedures. 
Statistical  evaluation  of  the  models  showed  higher  quality  of  mod¬ 
els  supported  with  variability  indices. 

The  cloud  cover  classifications  have  no  direct  connection  to  vari¬ 
ation  of  solar  irradiation  levels.  The  researchers  of  solar  energy 
systems  consider  the  solar  irradiation  as  an  input  signal  for  the 
state  of  system  operation.  Particular  days  are  grouped  on  the  basis 
of  signal  level  and  its  variation  in  time.  The  methods  for  classifi¬ 
cation  of  days  in  regard  to  solar  conditions  are  based  on  statistic 
instruments  and  give  different  results  depending  on  climatic  zones 
and  particular  purpose  of  classification.  The  classic  method  used 
for  classification  of  days  according  to  cloud  cover  is  based  on  the 
analysis  of  clearness  index  histograms  I<T. 

For  example,  Maafi  and  Harrouni  [16],  who  classified  solar  con¬ 
ditions  in  southern  Algeria,  used  two  threshold  criteria  in  their 
algorithm.  The  first  one  refers  to  frequency  distribution  of  global 
irradiance  (D),  and  the  second  one  to  the  beam  clearness  index). 
Three  classes  of  days  were  introduced: 


These  correlations  are  needed  in  calculations  of  radiation  sums  on 
the  inclined  surfaces,  when  only  total  irradiation  measurement  is 
available  for  the  horizontal  surface  [9].  Improvement  of  models  of 
Liu-Jordan  type,  from  one  hand  consists  in  testing  of  correlation 
in  shorter  and  shorter  time  intervals  [10,11],  and  from  the  other 
hand  in  introduction  of  additional  variables  like  the  sun  height, 
optical  mass,  atmosphere  humidity,  etc.  [12,13].  At  the  same  time, 
it  was  observed  that  under  the  cloudless  and  variable  cloudiness 
conditions,  the  daily  and  also  hourly  clearness  indices  often  take 
the  same  values.  It  results  then  in  poor  correlation  between  diffuse 
ratio  and  the  I<T  value.  A  substantial  improvement  in  model  quality 
can  be  obtained  by  introduction  of  clearness  variability  index  in 
hour  intervals  as  a  new  parameter  [14]. 

Parameters  of  solar  irradiation  variability  are  particularly  useful 
if  a  model  of  Liu-Jordan  type  is  to  be  based  only  on  measurements 
of  total  irradiation  on  the  horizontal  surface.  The  scientists  develop¬ 
ing  in  this  way  a  correlation  model  with  the  use  of  meteorological 


o  Clear  sky  day:  1  <  D  <  and  I<T  >  0.5. 
o  Partially  cloudy  sky:  Di  <  D  <  D2  and  I<T  >  0.5. 
o  Completely  cloudy  sky:  D2  <  D  or  D  <  D2  and  I<T  <  0.5. 

A  cumulative  distribution  function  was  used  to  determine  the 
thresholds  and  D2.  In  this  case  the  classification  was  carried  out 
for  the  need  of  photovoltaic  power  plants. 

Soubdhan  et  al.  [17],  who  investigated  climatic  conditions  in 
the  Caribbean,  introduced  another  method  for  classification  of  days 
according  to  daily  irradiation  conditions.  The  determination  of 
characteristic  types  of  daily  variation  was  related  to  the  dynamics 
of  energy  transformation  in  photovoltaic  systems.  The  basic  infor¬ 
mation  obtained  from  the  analysis  of  annual  results  was  the  daily 
frequency  histogram  of  the  beam  clearness  index  I<T  in  20  fractional 
ranges.  The  identification  of  histogram  classes  and  affiliation  prob¬ 
ability  of  particular  daily  histograms  for  the  established  class  was 
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performed  by  statistic  methods  based  on  a  finite  mixture  of  Dirich- 
let  distributions.  As  a  result,  the  authors  determined  three  features 
of  daily  irradiation: 

•  sunny  (S):  S  =  1  -  clear  sky  condition,  S  =  0  -  other, 

•  overcast  level  (C):  C  =  0  -  no  cloud;  C  =  1 ,  C  =  2  -  total  cloud  cover, 

•  dynamic  level  (D):  D  =  0,  D  =  l,  D  =  2  -  relates  to  the  velocity  of 
clouds  and  their  size. 

The  analysis  of  these  feature  combinations  led  to  the  identifica¬ 
tion  of  four  classes  for  daily  irradiation: 

(a)  First  class  -  for  S  =  l,  C  =  0  and  D  =  0.  This  class  is  represented 
by  a  clear  sky  with  low  number  of  clouds.  The  most  frequent 
are  periods  of  I<T  =  0.75.  In  the  climate  of  Guadalupe  7%  of  days 
were  affiliated. 

(b)  Second  class  -forS  =  l,C  =  l,D  =  l.  This  class  is  represented  by 
days  of  high  solar  irradiance  but  interrupted  with  medium  size 
clouds  of  limited  velocity. 

(c)  Third  class  -  for  S  =  0,  C  =  2,  D  =  0.  This  class  is  represented  by 
cloudy  days,  i.e.,  the  sky  is  covered  with  big  multilayer  clouds 
of  low  velocity.  The  frequency  of  Kj  is  highest  when  it  ranges 
between  0.1  and  0.2. 

(d)  Fourth  class  -  for  S  =  1 ,  C  =  2,  D  =  2.  This  class  is  represented  by 
days  of  sun  only  periodically  covered  with  small  clouds.  The 
index  frequency  I<T  reaches  two  peaks  -  for  0.25  and  0.75.  This 
class  predominates  in  this  climate  for  52%  days  in  a  year. 

The  classification  analyses  of  daily  irradiance  waveform  are  car¬ 
ried  out  to  distinguish  days  of  somehow  similar  features  within  the 
whole  population  of  days.  The  methods  and  classification  criteria 
assume  a  limited  number  of  classes;  as  a  result,  days  of  different 
irradiance  variation  are  included  in  the  same  class. 

A  method  for  classification  of  days  according  to  short-term 
variability  in  solar  irradiation  under  conditions  of  North  Europe  cli¬ 
mate  was  presented  by  the  Estonian  scientists  [18].  The  database 
included  a  period  from  April  to  September  of  1999-2002  (detailed 
analyses  concerned  the  May-July  period).  Ten  hours  that  were  sym¬ 
metrically  situated  around  the  solar  noon  were  considered.  The 
measuring  sequences  subjected  to  analysis  consisted  of  1-min  val¬ 
ues  of  total  radiation  G.  A  10-min  interval  was  taken  as  the  basic 
interval  for  investigation  on  irradiation  stability.  The  following 
index  of  radiation  intensity  variability  for  60  subsequent  intervals 
was  proposed: 


10 


If  this  index  value  exceeded  1000  W/(m2  min),  the  interval  was 
qualified  as  unstable.  It  was  assumed  that  if  the  unstable  inter¬ 
val  sum  exceeded  5h  (half  of  the  day),  this  day  was  arbitrarily 
counted  to  a  highly  unstable  group.  However,  if  the  unstable  inter¬ 
val  sum  was  equal  to  less  than  1  h,  this  day  was  conventionally 
regarded  as  stable  one.  In  summer  months  under  Estonia  climate 
the  cloudless  stable  days  amounted  to  10.7%  of  the  entire  ana¬ 
lyzed  period  and  they  supplied  1 6.7%  of  total  energy.  Days  of  strong 
cloud  cover,  that  were  also  qualified  as  stable  ones,  made  8.8%  of 
the  period  and  supplied  2.7%  of  total  energy.  At  the  next  stage  of 
investigations  there  was  developed  methodology  for  evaluation 
of  irradiation  conditions  variability  that  eliminated  the  notion  of 
basic  N-minute  interval.  Differences  AG  between  radiation  inten¬ 
sity  values  in  the  subsequent  minutes  of  the  day  were  presented  in 
the  form  of  distribution  curve.  The  curve  shape  was  characteristic 
for  variability  of  solar  irradiation  conditions.  Periods  of  gradients 
AG<  50  W/(m2  min)  predominate  on  the  cloudless  days.  On  the 


average,  16.2%  of  such  days’  length  were  characterized  by  AG  gra¬ 
dient  equal  to  zero.  However,  days  of  variable  cloud  cover  were  not 
univocally  parameterized. 

For  research  purposes,  the  irradiance  waveform  is  related  to 
1  h  or  sometimes  to  even  shorter  period,  e.g.  to  1  min.  It  can  be 
found  in  investigations  on  the  effect  of  short-term  irradiance  and 
its  components  in  classical  methods  of  calculations  of  global  radi¬ 
ation  on  differently  exposed  surfaces.  Short-term  solar  conditions 
vary  mainly  as  a  result  of  cloud  movement  and  that  is  why  dif¬ 
fuse  component  analysis  is  so  important.  It  appears  that  the  use  of 
short-term  data,  compared  to  hourly  averages,  does  not  affect  the 
results  of  monthly  global  radiation.  However,  if  a  useful  solar  frac¬ 
tion  is  determined,  i.e.  the  one  exceeding  the  established  threshold 
(utilizability  concept),  the  range  of  averaging  basic  measurement 
data  greatly  influences  the  shape  of  determined  waveforms.  The 
utilizability  based  on  3-min  averaging  is  by  10-30%  higher  for  crit¬ 
ical  values  ranging  from  200  to  500  W/m2,  when  compared  to  1-h 
averaging  [19].  Within  the  remaining  critical  value  ranges  these 
differences  are  not  so  big.  Simulations  of  solar  energy  conversion 
efficiency  in  power  systems  based  on  short-term  radiation  data 
should  also  lead  to  correction  of  the  results  based  on  hourly  aver¬ 
ages. 


3.  Methodology 


In  simulation  computations  for  solar  system  heating  efficiency 
and  testing  of  real  systems  the  hourly  average  irradiation  is  used  as 
a  standard.  The  comparison  of  efficiency  and  dynamics  of  heating 
solar  systems  within  the  energy  diagnostics  requires  definite  evalu¬ 
ation  of  radiation  conditions  during  tests.  It  is  advised  to  introduce 
a  parameter  (an  indicator)  to  identify  the  days  of  different  cloud 
covers  in  the  same  group  of  daily  global  irradiation.  The  effects  of 
heating  solar  system  operation  depend  not  only  on  daily  dose  but 
also  on  its  waveform  during  a  day.  The  authors  introduce  a  concept 
of  daily  parameterization  for  solar  irradiance  and  present  two  dif¬ 
ferent  methods  to  determine  the  variation  parameter  without  any 
preference.  Both  of  them  have  their  advantages  and  disadvantages, 
and  in  each  of  the  methods  a  particular  day  is  definitely  described  by 
means  of  rational  indexes,  however,  sometimes  of  different  values. 

The  first  method  uses  an  elementary  index  of  property  varia¬ 
tion  as  a  standard  deviation  a  related  to  the  mean  value  (variation 
indicator  Wz).  The  investigated  property  is  global  solar  irradiance  G 
recorded  at  minute  intervals.  Solar  days  in  this  method  are  divided 
into  equal  time  intervals  Ar  of  k-minutes.  The  number  of  intervals 
n  results  from  the  number  of  m  minutes  from  sunrise  to  sunset: 


Range  variation  indices  are  determined  for  particular  day: 

(V\4)AT=aATGAT  (2) 

where  AGAt  is  average  irradiance  in  a  time  interval,  W/m2. 

Variation  index  for  the  whole  day  ( Wz)d  is  weighted  average  of 
time  intervals’  values  related  to  arithmetic  mean  irradiance  for  the 
whole  day.  The  weight  is  the  share  of  irradiation  from  i-interval  in 
the  daily  sum  H  that  is  expressed  in  this  paper  in  Wh/m2: 


(Wz)d  = 


(ELi^Ca^)  /60 H 
60  H/m 


Introduction  of  weight  for  intervals’  values  aimed  at  intensifi¬ 
cation  of  variation  rank  of  radiation  conditions  in  the  periods  of 
highest  radiation  intensity.  This  method  is  thus  useful  for  precise 
determination  of  solar  conditions  for  heating  systems,  where  the 
periods  of  low  radiation  are  not  important.  The  value  of  ( Wz)d  is 
strongly  influenced  by  the  accepted  daily  number  of  intervals.  The 
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Fig.  1.  Selected  solar  irradiance  waveform  composed  of  the  constant  component  and  subsequent  harmonics  on  a  background  of  real  waveform. 


interval  duration  k  should  depend  on  the  dynamics  of  conversion 
phenomena  in  a  given  type  of  objects.  The  criterion  can  be,  e.g. 
thermal  time  constant.  It  is  determined  for  collectors  according  to 
standard  procedures  (EN  12975  [20]  and  ISO  9806  [21]).  Obtained 
thermal  time  constants  depend  on  collector  structure  and  range 
from  few  minutes  for  flat  plate  collectors  to  several  minutes  for 
heat  pipe  tube  collector  sets.  Time  constants  for  heating  systems, 
where  the  collectors  cooperate  with  heat  exchangers  and  volume 
tanks  of  high  thermal  inertia,  amount  to  tens  of  minutes.  This  paper 
presents  daily  parameterization  of  solar  irradiation  with  the  use  of 
( Wz)d  index  for  the  most  frequently  used  meteorological  data  range 
of  60  min,  with  introduction  of  ( Wz)d60. 

The  second  evaluation  method  treats  solar  radiation  intensity  as 
a  time  varied  signal  h{  r).  Analogically  to  electric  signals,  the  signal 
waveform  has  been  analyzed  with  the  use  of  Fast  Fourier  Transform 
(FFT).  Harmonic  components  ft;  and  their  carried  power  have  been 
determined.  The  daily  radiation  spectrum  analysis,  carried  out  for 
several  years,  shows  that  main  signal  power  is  carried  by  the  con¬ 
stant  component  h0  (mean  value)  and  the  two  first  harmonics  h\ 
(of  24-h  period)  and  h2  (of  12-h  period)  [22]. 

Fig.  1  shows  the  comparison  between  waveforms  composed  of 
the  constant  and  harmonic  components  for  a  particular  day.  It  is 
clear  that  application  of  harmonic  component  higher  than  the  sec¬ 
ond  one  causes  the  signal  deformation.  The  real  sun  trajectory  in 
relation  to  measurement  surface  is  represented  by  the  waveform 
consisted  of  the  constant  component  and  the  two  first  harmonics. 
The  higher  harmonics  can  be  regarded  as  interference  in  the  signal 
waveform  spectrum  h{ r). 

That  is  why,  analogically  as  in  electronics,  the  index  THD  ( Total 
Harmonic  Distortion)  has  been  established  for  the  share  of  harmon¬ 
ics  that  are  qualified  as  interference  in  the  full  signal  spectrum: 

_ _ V^32  +  ^4 2  +  ^52  H  F  hi7  ^ 

\f ho  +  h\  +  h2  +  h2  ~-\-h 4  +  •  •  •  +  hfc 

where  h0  constant  component  value,  and  hi  amplitudes  of  subse¬ 
quent  harmonics. 

It  follows  from  frequency  analysis  of  the  h{  r)  signal  that  ampli¬ 
tude  values  of  particular  harmonics  depend  on  cloud  cover  levels, 
and  their  percentage  share  in  the  spectrum  depends  on  frequency 
of  cloud  cover  changes.  At  higher  cloud  cover  the  number  of 


Fig.  2.  Deterministic  solar  irradiance  waveform  of  the  tested  surface. 

harmonic  components  increases,  but  the  values  of  main  component 
amplitudes  decrease. 

The  method  for  daily  solar  irradiance  parameterization  by 
means  of  THD  index  does  not  directly  take  into  account  the  varia¬ 
tion  weight  in  the  hours  of  highest  doses,  contrary  to  the  method 
with  ( Wz)d  index. 


Time 


Fig.  3.  Solar  irradiance  variation  on  a  day  of  full  cloud  cover. 
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4.  Exemplary  results 

The  irradiance  waveform  database  used  for  testing  of  the 
described  parameterization  methods  includes  several  hundred 
days.  The  measurement  surface  is  south  oriented  and  tilted  at  35°. 
During  sunny  days  (clear  sky)  the  variation  index  of  solar  irradi¬ 
ance  results  from  deterministic  position  of  the  sun.  Fig.  2  shows 
solar  irradiance  waveform  (for  official  time)  during  a  clear  sky  day 
of  high  atmosphere  visibility.  The  variation  indices  are  then  equal 
to  ( Wz)d60  =  0.046,  and  THD  =  0.050. 

Under  conditions  of  almost  full  cloud  cover,  when  periodical 
global  radiation  results  from  diffuse  radiation,  irradiance  variation 
is  also  observed  (Fig.  3).  The  absolute  variation  range  is  small  then, 
but  relative  variations  are  substantial  and  thus  the  indexes  THD  and 
(Wz)d  60  reach  values  several  times  higher  than  that  at  deterministic 
radiation  variation.  The  evaluation  of  irradiance  waveforms  during 
days  of  continuous  overcast  has  only  methodical  meaning  because 
solar  heating  systems  are  not  active  then. 

The  subsequent  graph  (Fig.  4)  shows  the  solar  irradiance  wave¬ 
form  on  a  day  of  several  cloud  cover  phases.  Until  11:15  the 
clouds  did  not  cover  the  sun,  but  irradiance  changed  quickly 
and  in  a  wide  range  in  two  subsequent  hours.  There  was 


Time 


Fig.  4.  Solar  irradiance  variation  on  a  day  of  varied  overcast. 


no  beam  radiation  between  13:30  and  15:15  on  the  mea¬ 
surement  surface.  Then,  until  the  end  of  the  day,  overcast 
varied.  The  parameters  then  were  computed  as  THD  =  0.5 10,  and 
(Wz)d60  =  0.476. 


* 


a> 


■a 

S3 

S3 

O 

a 

> 


u 

'o 

C/5 


Number  of  sample 

Fig.  5.  The  comparison  of  daily  indices  of  solar  irradiance  variations:  THD  and  (Wz)d60. 
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Fig.  6.  Solar  irradiance  conditions  parameterized  differently  by  indexes  THD  and  ( Wz)d6 o- 
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Fig.  7.  The  variation  of  test  surface  irradiance  on  19.08.2002. 
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Fig.  9.  The  variation  of  test  surface  irradiance  on  02.07.2005. 


Fig.  8.  The  variation  of  test  surface  irradiance  on  01.08.2005. 


The  best  convergence  of  proposed  indexes  occurs  in  the  case 
when  for  (Wz)d,  interval  Ar  =  60min.  Fig.  5  presents  the  compari¬ 
son  of  both  indices  for  over  1 000  days  of  the  period  between  80  days 
before  and  80  days  after  the  summer  solstice  (i.e.  between  1  st  Apr. 
and  1 0th  Sept.).  During  the  remaining  periods  of  year  the  test  stand 
is  shaded  and,  thus,  the  comparison  results  are  not  reliable.  It  has 
been  found  that  at  THD  values  less  than  0.3,  the  index  ( Wz)^ 60  is  on 


the  average  by  18.6%  lower  in  the  relative  scale,  at  difference  stan¬ 
dard  deviation  of  14.1%.  Within  the  THD  range  between  0.3  and  0.7, 
the  index  (Wz)d60  is  on  the  average  higher  by  0.3%,  but  difference 
standard  deviation  amounts  to  20.8%. 

In  some  special  conditions  of  solar  irradiance  waveform,  distinct 
differences  can  be  observed  on  the  test  surface.  Such  example  is  a 
generally  clear  day,  when  radiation  was  limited  by  fog  during  morn¬ 
ing  and  evening  hours  (Fig.  6).  The  hourly  indexes  ( Wz)h  changed 
in  the  wide  range  from  1 .1 28  at  9  o’clock  to  0.006  at  13  o’clock.  The 
arithmetic  mean  of  hourly  indexes  amounted  to  0.3,  but  the  daily 
index  (Wz)d60,  calculated  with  the  use  of  particular  hours’  rank, 
was  reduced  to  0.066.  The  irradiance  waveform  variation  parame¬ 
terized  by  THD  index  was  equal  to  0.245.  That  was  the  case,  when 
both  indexes  differed  distinctly.  In  such  unusual  atmospheric  con¬ 
ditions  the  index  (Wz)d60  better  characterized  irradiance  variation 
in  the  context  of  solar  energy  use  in  thermal  processes. 

The  subsequent  graphs  (Figs.  7-9)  indicate  that  basing  on  daily 
global  irradiation  one  cannot  conclude  on  irradiance  waveforms 
during  a  day.  Three  exemplary  days  have  been  selected  from 
the  database  (19th  Aug.,  1st  Aug.,  2nd  Jul.),  of  almost  identical 
global  sum  equal  to:  6328  Wh/m2,  6366  Wh/m2  and  6267  Wh/m2, 
respectively.  However,  variation  parameterization  showed  distinct 
differences  for  those  days:  on  the  first  (clear)  day  (Wz)d60  =  0.069 
and  THD  =  0.086,  on  the  second  day  0.307  and  0.317,  and  on  the 
third  day  0.363  and  0.353,  respectively. 


Fig.  10.  The  waveform  of  tank  heating  with  battery  of  collectors  on  a  background  of  solar  irradiance  on  tested  surface  on  02.07.2005. 
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The  substantial  difference  in  solar  radiation  waveforms  on 
the  compared  days  (Figs.  7-9)  resulted  in  exploitation  effects  of 
exemplary  solar  heating  system,  consisted  of  tubular  heat  pipe 
collectors  and  a  volume  heat  exchanger  [23].  These  days  pro¬ 
vided  similar  water  initial  temperature  in  the  accumulation  tank, 
while  water  disposal  during  solar  system  activity  was  small  and 
almost  identical.  On  2nd  Jul.  between  10:00  and  16:00,  11  pauses 
of  medium  flow  were  observed  (lasting  together  22  min),  while 
the  circulation  pump  medium  flow  amounted  to  2085  dm3.  On  1st 
Aug.  between  10:00  and  16:00,  12  pauses  were  observed  (lasting 
together  16  min),  while  the  circulation  pump  medium  flow  was 
equal  to  2140  dm3.  On  19th  Aug.  the  circulation  pump  was  work¬ 
ing  continuously  and  its  flow  was  equal  to  2315  dm3.  As  a  result, 
the  daily  efficiency  of  collectors  amounted  to  28.6%  on  the  first 
day,  31.5%  on  the  second  day,  and  38.5%  on  the  third  day.  More¬ 
over,  it  should  be  mentioned  that  vacuum  collector  section  of  the 
investigated  system  works  at  high  temperature  in  the  accumulation 
tank. 

The  variability  of  solar  irradiance,  particularly  in  the  middle  part 
of  the  day,  influences  the  waveform  of  tank  heating  by  the  set  of  col¬ 
lectors.  Fig.  10  presents  the  process  of  medium  heating  up  in  the 
battery  of  vacuum  collectors  and  of  hot  water  in  the  tank  on  the 
day  characterized  by  high  variation  indices.  Solar  irradiance  varied 
very  much  between  1 1 :20  and  13:05  and  then  between  14:20  and 
14:50  as  well  as  from  15:30  to  1 6:00.  In  the  first  interference  period, 
the  dynamics  of  heating  slowed  down  and  stopped  in  the  second 
period.  The  highest  temperature  increase  occurred  between  13:05 
and  14:20,  i.e.  in  the  period  of  directly  visible  sun’s  disk  and  the 
increased  radiation  intensity  due  to  reflection  from  clouds.  How¬ 
ever,  the  final  solar  irradiation  conversion  efficiency  computed  for 
the  whole  day  was  low. 

5.  Conclusions 

In  simulations  of  solar  heating  and  photovoltaic  systems  the 
standard  averages  of  hourly  radiation  data  are  used.  The  more 
professional  simulations  require  short-term  data  (even  1-min 
data);  this  is  particularly  important  in  control  system  problems. 
Operational  tests  of  real  objects  aimed  at  verification  of  design 
calculations  and  comparisons  between  solar  system  solutions.  The 
successful  completion  of  these  tests  calls  for  precise  determination 
of  experimental  conditions.  The  global  irradiation  data  are  not  pre¬ 
cise  enough  and  without  additional  parameterization  of  variation 
are  of  a  limited  practical  value. 

The  parameterization  of  daily  solar  irradiance  variation  can 
be  performed  with  the  use  of  classic  indexes:  mathematical  THD 
and  statistical  one  (based  on  the  concept  of  standard  deviation). 
The  concepts  of  both  indices  have  different  methodological  back¬ 
grounds  but  provide  similar  results,  when  a  basic  range  in  the 
statistical  method  amounts  to  60  min.  Clear  days  of  solar  irradiance 
deterministic  variation  are  parameterized  on  the  tested  surface  by 
indexes  of  the  following  values:  THD  from  0.05  to  0.08  and  ( Wz)d6o 
between  0.05  and  0.07.  The  increasing  index  values  point  out  at 


higher  and  higher  share  of  stochastic  interference.  The  advantage 
of  THD  index  is  its  standardization  that  limits  its  values  within  the 
(0.1 )  range.  The  ( Wz)d60  index  value  can  be  higher  than  1 ,  but  it  takes 
into  account  higher  irradiance  values  in  the  better  way  (through  the 
introduction  of  weight  parameter),  which  is  particularly  advanta¬ 
geous  for  solar  heating  systems. 

Days  of  similar  global  irradiation  can  be  characterized  by  dif¬ 
ferent  variations  of  irradiance  in  subsequent  hours  and  that  is  why 
the  introduction  of  THD  or  (Wz)d60  provides  good  effects  for  the 
purpose  of  identification.  The  problem  of  selection  between  the 
proposed  indices  requires  further  research. 
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